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Abstract 

The acceleration of interstellar pick-up ions as well as solar wind species has been 
observed at a multitude of interplanetary (IP) shocks by different spacecraft. The ef- 
ficiency of injection of the pick-up ion component differs from that of the solar wind, 
and is expected to be strongly enhanced at highly oblique and quasi-perpendicular 
shock events, in accord with inferences from in situ observations. This paper explores 
theoretical modeling of the phase space distributions of accelerated ions obtained 
by the Ulysses mission for the Day 292, 1991 shock associated with a corotating 
interaction region, encountered before Ulysses' fly-by of Jupiter. A Monte Carlo 
simulation is used to model the acceleration process, adapting a technique that 
has been successfully tested on earlier IP shocks possessing minimal pick-up ion 
presence. Phase space distributions from the simulation technique for various low 
mass ions are compared with SWICS and HI-SCALE data to deduce values of a 
"turbulence parameter" that controls the efficiency of injection, and the degree of 
cross-field diffusion. Acceptable fits are obtained for the H + and He + populations 
using standard prescriptions for the pick-up ion distribution; He ++ spectral data 
was only fit well for scenarios very close to the Bohm diffusion limit. It is also found 
that the simulation successfully accounts for the observation of energetic protons 
farther upstream of the forward shock than lower energy pick-up protons, using the 
same turbulence parameter that is required to achieve reasonable spectral fits. 

Key words: Interplanetary shocks, Co-rotating interaction regions, Ulysses 
mission, Diffusive shock acceleration, Shock turbulence 
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1 Introduction 



Particle acceleration at collisionless shocks is believed to be a common 
phenomenon in space plasmas in a diversity of environments, ranging from 
the inner heliosphere to the central regions of distant galaxies. In the helio- 
sphere, evidence to support the belief that such a mechanism can efficiently 
produce non-thermal particles includes direct measurements of accelerated 
populations in various energy ranges at the Earth's bow shock (e.g. Scholer 
et al. 1980; Mobius et al. 1987; Gosling et al. 1989) and interplanetary shocks 
(for the pre-Ulysses era see, for example, Sarris and Van Allen 1974; Gosling 
et al. 1981; Decker, Pesses and Krimigis, 1981; Tan et al. 1988). The motiva- 
tions for developing theories of shock acceleration are therefore obvious, and 
a variety of approaches have emerged. One possible means for the generation 
of non-thermal particles is the Fermi mechanism, often called diffusive shock 
acceleration; it is this process that is the focus of this paper. 

There are various approaches to modelling diffusive shock acceleration. 
Among these are convection-diffusion differential equation approaches (e.g. 
Kang and Jones 1995), and the kinematic Monte Carlo technique of Ellison 
and Jones (e.g., see Jones and Ellison, 1991, and references therein), which de- 
scribes the injection and acceleration of particles from thermal energies; both 
of these can address spectral and hydrodynamic properties. This latter simu- 
lational approach is the central tool for the analysis here, and is ideally suited 
to the interpretation of the modest time-resolution shock data acquired by 
Ulysses. Hybrid and full plasma codes (e.g. Quest 1988; Burgess 1989, Winske 
et al. 1990; Trattner and Scholer 1991; Giacalone, Burgess and Schwartz 1992; 
Liewer, Goldstein, and Omidi 1993; Kucharek and Scholer 1995) provide con- 
trasting probes of shock environs, with an emphasis primarily on plasma struc- 
ture and wave properties in the environs of shocks. 

Comparisons of predictions from theoretical models with observed phase- 
space distribution data are informative probes. The first detailed theory/data 
comparison along these lines was performed by Ellison, Mobius and Paschmann 
(1990) in the case of the quasi-parallel portion of the Earth's bow shock, com- 
paring predictions of the Monte Carlo method with particle distributions of 
protons, He ++ and a C, N and O ion mix obtained by the AMPTE exper- 
iment. The agreement between model predictions and data was impressive. 
Ellison, Mobius and Paschmann (1990) concluded that a successful fit was 
possible only in the non-linear acceleration regime, when the dynamic effects 
of the accelerated particles are crucial to the determination of the shock struc- 
ture. Since this pioneering work, successful comparisons of other theoretical 
techniques with data from the Earth's bow shock have been performed. These 
include the hybrid plasma simulations of Trattner and Scholer (1991) and Gi- 
acalone et al. (1993), and solutions to the convection-diffusion equation (Kang 
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and Jones 1995), both of which have yielded good agreement with both the 
bow shock data and the Monte Carlo technique (Ellison et al. 1993). 

Such a comparison between theory and experiment was extended to the 
domain of interplanetary shocks in the work of Baring, Ogilvie, Ellison and 
Forsyth (1997, hereafter BOEF97; see also Kang & Jones 1997, for application 
of their convection-diffusion equation technique). In this development, impres- 
sive agreement was found between the Monte Carlo predictions and spectral 
data obtained by the Solar Wind Ion Composition Spectrometer (SWICS) 
aboard Ulysses, in the case of two shocks observed early in the Ulysses mis- 
sion. Such agreement was possible only with the assumption of strong particle 
scattering (i.e. near the Bohm diffusion limit) in the highly oblique candidate 
shocks. For a third shock, detected a month later, the comparison failed, with 
significant differences arising in the 500-800 km/sec range of the phase space 
distribution. BOEF97 attributed this discrepancy to the omission of pick-up 
ions from the model: such an extra component would be expected to provide a 
substantial contribution to the accelerated population in this particular event. 

This paper addresses the role of pick-up ions in such shocks via modeling 
the accelerated population for the specific interplanetary shock observed on 
day 292 of 1991, detected by the SWICS and HI-SCALE instruments aboard 
the Ulysses spacecraft at around 4.5 AU, as reported in Gloeckler et al. (1994). 
Phase space distributions from the simulations are compared with SWICS and 
HI-SCALE data, yielding acceptable fits for the H + and He ++ populations, 
using standard prescriptions for the injected pick-up ion distribution, by ad- 
justing a single turbulence parameter rj . Using this same rj , the simulation 
results also successfully account for the observation of energetic protons far- 
ther upstream of the forward shock than lower energy pick-up protons, since 
a rigidity-dependent diffusion is used in the models. 



2 The Monte Carlo Simulation Technique 

The Monte Carlo simulation technique used here has been invoked in pre- 
vious applications to shocks in the heliosphere, as discussed above, and is 
extensively described in a number of papers (Ellison, Jones, & Eichler 1981; 
Ellison, Jones & Reynolds 1990; Jones & Ellison 1991; Baring, Ellison & Jones 
1994; Ellison, Baring, & Jones 1996). It follows closely Bell's (1978) test par- 
ticle approach to diffusive acceleration. Particles are injected upstream and 
allowed to convect into the shock, meanwhile diffusing in space so as to effect 
multiple shock crossings, and thereby gain energy through the shock drift and 
Fermi processes. For the interplanetary shock that is the focus of this paper, 
throughout u\ and u<i shall denote the upstream and downstream flow speed 
in the shock rest frame, respectively. 
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The particles gyrate in a laminar electromagnetic field. Their trajectories 
are obtained by solving a fully relativistic Lorentz force equation in the shock 
rest frame, in which there is, in general, a u x B electric field in addition to 
the magnetic field. The effects of magnetic turbulence are modeled by scat- 
tering these ions in the rest frame of the local fluid flow. These collisions, 
effectively mimicking diffusive transport instigated by Alfven waves, are as- 
sumed to be elastic in the local fluid frame, an assumption that is valid so 
long as the flow speed far exceeds the Alfven velocity. Otherwise, the scatter- 
ing centers are not anchored in the fluid frame, and the collisions are inelastic, 
thereby introducing significant second-order Fermi (stochastic) acceleration. 
For interplanetary shocks in general, and in particular for the Day 292, 1991 
shock considered here, the elastic scattering assumption is somewhat violated 
due to a low shock Alfvenic Mach number. Yet the stochastic acceleration 
contribution is expected to be relatively small, perhaps at most of the order 
of 25%, so its inclusion is deferred to later work. 

The simulation can routinely model either large-angle or small-angle scat- 
tering. In this paper, large-angle scattering is employed, motivated by previous 
observations (e.g., Hoppe et al. 1981; Balogh et al. 1993) indicating very tur- 
bulent fields in IP shocks. Between scatterings, particles are allowed to travel 
a distance in their local fluid frame that is exponentially distributed about the 
scattering mean free path A , with: 



Here r g = pcj (qB) is the gyroradius of an ion of momentum p = mv , mass 
m , and charge q . Also r g \ = mu\, x cj {qB) is the gyroradius of an ion with a 
speed v equal to the far upstream flow speed normal to the shock plane, u\, x , 
where x denotes the direction normal to the shock plane. The mean free path 
scale Ao is set proportional to r g i with constant of proportionality r\ such 
that A = i]r gl . Following previous Monte Carlo work, for simplicity we set 
a — 1 , a specialization that is appropriate for interplanetary plasma shocks 
(e.g., see Ellison et al. 1990; Mason, Gloeckler & Hovestadt 1983; Giacalone, 
Burgess & Schwartz 1992 for discussions on the micro-physical expectations 
for a ) , yet is easily generalizable in the simulation. Since A > r g for physically 
meaningful diffusion resulting from gyro-resonant wave-particle interactions, 
the a = 1 case is also motivated on fundamental grounds. 

Cross-field diffusion emerges naturally from the simulation, since at every 
scattering, the direction of the particle's momentum vector is randomized in 
the local fluid frame, with the resulting effect that the gyrocenter of a par- 
ticle is shifted randomly by a distance of the order of one gyroradius in the 
plane orthogonal to the local field. Transport perpendicular to the field is 
then governed by a kinetic theory description, so that the ratio of the spa- 
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tial diffusion coefficients parallel (k\\ — Xv/3 ) and perpendicular (k±) to the 
mean magnetic field is given by k±/ k\\ = 1/(1 + 1] 2 ) (see Forman, Jokipii & 
Owens 1974; Ellison, Baring & Jones 1995, for detailed expositions). Hence, 
i] couples directly to the amount of cross-field diffusion, and is a measure not 
only of the frequency of collisions between particles and waves, but also of the 
level of turbulence present in the system, i.e. is an indicator of (5B/B) . Note 
that 1] = 1 is the Bohm diffusion limit of quasi-isotropic diffusion, presumably 
corresponding to (5B/B) ~ 1 . As will become evident in the Section 3, by 
virtue of its connection to cross-field diffusion, 77 plays an important role in 
determining the injection efficiency of low energy particles. We note that other 
cross-field transport effects such as that incurred by so-called field-line wan- 
dering (e.g. see Giacalone & Jokipii 1999) can be incorporated parametrically 
in the Monte Carlo technique, and will be addressed in future work. 

While i] and a serve as parameters here, in principal they can be cal- 
culated for a given wave field. However, in practice, precision is limited by 
in situ magnetometer data, and exact determination of 3-dimensional parti- 
cle diffusion properties from a 1-dimensional field data stream is impossible. 
Yet insights can be gained by using established theoretical formalism, such 
as random fluctuation theory (e.g. Jokipii & Coleman 1968), or equivalently 
quasi-linear diffusion theory (see, e.g., Kulsrud & Pierce 1969; Jones, Birming- 
ham, & Kaiser 1978), to estimate components of the diffusion tensor along and 
orthogonal to the mean field. These components can be expressed in terms of 
integrals of components of the power spectral tensor for field fluctuations. Such 
slab models for turbulence can be extended to consider higher-dimensional 
field fluctuations (e.g. see Bieber et al. 1994). Guidance for the behavior and 
values of a and i] can be gleaned from such investigations, and an analysis 
along these lines to facilitate the interpretation of Ulysses IP shock spectral 
data will be the subject of future work. 

The simulation technique makes no distinction between accelerated parti- 
cles and thermal ones. Therefore the injection efficiency is determined solely 
by the a particle's ability to diffuse back upstream after it has encountered 
the shock for the first time. In the case of the Earth's bow shock (Ellison et 
al. 1990), this permits the development of a nonlinear model that includes the 
effects of the accelerated particles on the dynamics of the shock. However, in 
the relatively weak IP shock considered here, a steep power-law is present, and 
the poor injection efficiency reduces the number density of energetic particles, 
eliminating nonlinear hydrodynamic effects. Upstream plasma quantities are 
input from observational data, as discussed below, and downstream quantities 
are determined using the full MHD Rankine-Hugoniot relations. The simu- 
lation is capable of producing energy spectra/fluxes, and hence phase space 
distributions, at any location upstream or downstream of the shock and in 
any reference frame including that of a spacecraft. This makes the simulation 
ideal for comparison with observational data. 
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3 Modeling the Ulysses Event of Day 292, 1991 

Our case study focuses on the forward shock of a CIR encountered by 
Ulysses on Day 292 of 1991, for which downstream particle distributions were 
obtained for comparison with data published in Gloeckler et al. (1994). To 
determine parameters for this shock that are appropriate input for the Monte 
Carlo simulation, the observational results detailed in Gloeckler et al. (1994) 
are used along with the data compilations of Balogh et al. (1995), Gonzalez- 
Esparza et al. (1996), and Hoang et al. (1995). These sources indicate that the 
sonic Mach number M s of the shock was 2.53, with an unstated uncertainty. 
This controls the temperature of the thermal solar wind protons, but does 
receive contributions from electrons and alpha particles. Other parameters 
include the angle 9 Bn i = 50° ±11° the upstream magnetic field makes with 
the shock normal, B x = 0.8 nT, and magnetic compression ratio B 2 /Bi = 
2.50±0.13 . Here B\ and B 2 are the total magnetic field magnitudes upstream 
and downstream of the shock, respectively. In addition, the normalization of 
solar wind distributions was established using n p = 2.0cm~ 3 as the solar 
wind proton density, and n He = 4.0 x 10~ 2 cm -3 as that for solar wind He ++ 
(within the uncertainties of values listed Table 1 in Gloeckler et al. 1994). 

Although fluctuations in the plasma data preclude precise determination 
of the velocity compression ratio of the shock, Gloeckler et al. inferred a value 
of r = u\/u2 = 2.4 ± 0.3. This inference is somewhat inconsistent with the 
observed t>~ 5 - 75 power-law behavior of the phase space density at high ener- 
gies, which suggests a compression ratio nearer 2. The choice of compression 
ratio is controlled by the shock Mach number M T = M S M A / ^ M| + M| , re- 
ceiving contributions from sonic and Alfvenic Mach numbers, M s and M A 
respectively, which can be written in the form (e.g. Baring et al. 1997): 

M s w 8.52 . Ul0 ° , M A « 4.58 ^ , (2) 

^T p4 + T e4 5 -5 V 

where wioo is the shock speed u± tX in units of 100 km/sec, T p4 and T e4 are 
the proton and electron temperatures in units of 10 4 K, £>„ 5 is the upstream 
field strength in units of 10~ 5 Gauss, and n p is the proton density in cm -3 . 
Using upstream values of T e4 = 2.59 and T p4 = 1.67, as tabulated in Hoang 
et al. (1995), we find that for U\ >x = 55 km/sec, M s 2.5. Adding in the 
magnetic field gives M T ps 2.1 . To yield a compression ratio r = (7 + l)/[y — 
1 + 2/M|] =2.1 for adiabatic index 7 = 5/3, one must have a slightly lower 
value of M T p^ 1.85 (i.e. perhaps higher field) . In this paper, we adopted a 
value of r = 2.1 , which lies within the margin of error of Gloeckler et al.'s 
inferred value, to generate power-law slopes commensurate with the data. Due 
to significant observational uncertainties, ranges of parameters are permissible. 



6 



In the simulation runs, the temperature T a of the a particles was cho- 
sen to be four times that of the solar wind protons. No guidance for this 
assumed upstream temperature was provided in Gloeckler et al. (1994), nor 
in the aforementioned Ulysses data compilation papers. However, in the IP 
shock theory/data comparison of Baring et al. (1997), for the one shock with 
good measurements of the He ++ distribution, namely the 91097 event, the 
upstream a particle thermal velocity was comparable to (actually slightly 
greater) than that of the solar wind protons, implying T a rs 4T p . Adopting 
the same temperature ratio T a /T p in this work was motivated by this prece- 
dent, and the expectation that the adiabatic expansion of the solar wind from 
2.7AU to 4.5AU should roughly preserve the ratios of the temperatures per 
nucleon for these two ions. Note that the correction to the sonic Mach number 
in Eq. (2) incurred by solar wind helium is then approximately 4%. 

To facilitate transformation of simulation results measured in the shock 
rest frame to the approximate spacecraft frame, so as to compare directly with 
Gloeckler et al.'s published distributions, a value for the solar wind plasma 
speed of v sw ^ = 393 ± 12 km/s was assumed after the shock encounter. This 
identifies the downstream fluid frame. No values for the shock speed were listed 
in Gloeckler et al. (1994), as is common practice in other expositions (e.g. see 
Burton et al. 1992; Baring et al. 1997). Hence selection of a value for u\ tX 
consistent with other data became necessary. From the long-term solar wind 
radial velocity traces in Figure 2 of Gonzalez-Esparza et al. (1996) we can 
set a lower bound on the solar wind speed of i> STOj i > 330 km/s prior to the 
shock passage, i.e. upstream. Assuming collinear velocities, i.e. that the solar 
wind is not deflected by the shock on large scales, we have the speed relations 
ui, x = v sh - v SWjl and u 2yX = v sh - v sw , 2 in addition to r = u hx /u 2)X , where 
v s h is the unknown shock speed in the spacecraft frame. These solve to yield 
Ui iX ~ 55 km/s and u 2<x ~ 28 km/s for r = 2.1 for v SWt i ~ 360 km/s; higher 
values of u\ >x and u 2 ^ x would be realized for smaller choices of v SW; i . 

The other input for the Monte Carlo simulation is that for the incoming 
pick-up ion distributions. Gloeckler et al. (1994) used a simplified form for 
these distributions; here we adopt the slightly more developed expression used 
in Ellison, Jones & Baring (1999) that is modeled on the seminal work of 
Vasyliunas & Siscoe (1976). This form incorporates the gravitational focusing 
of interstellar neutrals, the physics of their ionization as a function of distance 
from the sun (approximately at solar minimum), and adiabatic losses incurred 
during propagation away from the sun, and is similar in conception to pick-up 
ion distributions used in le Roux, Potgieter & Ptuskin (1996). Our pick-up ion 
model therefore provides both the detailed shape and normalization of these 
superthermal distributions, which are distinctly different for H + and He + , 
largely due to the patently different ionization and charge exchange rates for 
their corresponding interstellar neutrals. 
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Ulysses IP shock on Day 292, 1991 




L°gio[ v /v sw ] 

Fig. 1. Phase space velocity distribution functions for data collected by the 
Ulysses mission for the shock on Day 292 of 1991, specifically for the interval 
1991.292.0400-293.0400. The velocity axis is the ratio of the ion speed v, as mea- 
sured in the spacecraft frame, to the solar wind speed. The data are for H + (filled 
circles for SWICS data; open circles for HI-SCALE points) solar wind and pickup 
ions, He + (filled triangles) pickup ions, and He ++ (filled squares), mostly so- 
lar wind ions, and are taken from Gloeckler et al. (1994). The histograms are the 
corresponding Monte Carlo models of acceleration of these species (heavy weight 
corresponds to solar wind ions) for u\ = 55 km/sec, using the optimal choice of 
plasma shock parameters from Gloeckler et al. (1994) and sources indicated in the 
text. The model assumed r] = X/r g = 5.5 and a shock of compression ratio r = 2.1 , 
corresponding to diffusive acceleration power-laws of index —5.73, is indicated. 

Downstream distributions for thermal and accelerated ions from the Monte 
Carlo simulation are compared with SWICS and HI-SCALE measurements 
taken in the frame of the spacecraft on the downstream side of the Day 292, 
1991 shock in Fig. 1. The Ulysses data are those exhibited in Fig. 1 of Gloeck- 
ler et al. (1994), and the simulation data are transformed to the spacecraft 
frame as described above. The solar wind and pick-up ion parameters and 
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abundances are fairly tightly specified, so that the theoretical model has one 
largely free parameter, the ratio of the particle mean free path to its gyrora- 
dius, i] = \/r g . As noted in Ellison, Baring & Jones (1995) and Baring et 
al. (1997), the efficiency of acceleration of thermal ions in oblique shocks, i.e. 
the normalization of the non-thermal power-law, is sensitive to the choice of 
i] , so this parameter was adjusted to obtain a reasonable "fit" to the Ulysses 
data. Ellison, Jones & Baring (1999) note that for their application to the 
solar wind termination shock, the efficiency is not as sensitive to i] when su- 
perthermal pick-up ions enter the problem, a property that we also replicate 
here. Here, for r\ > 3 , the accelerated pick-up ion phase space density is orders 
of magnitude above that of the thermal ions, for both H and He . 

In Fig. 1, the theory/data comparison is overall not quite as good as the 
ones in Baring et al. (1997) for Ulysses IP shock data at around 2-3 AU, where 
pick-up ions are less of a factor. Yet the fits here do model the accelerated 
pick-up ions very well, for rj = 5.5 , a value that is slightly higher than those 
inferred in the fits of Baring et al. (1997), but is still consistent with a moderate 
level of field turbulence (i.e. here k±/k\\ ~ 0.03). Note that since the inferred 
value of i] = 5.5 for H + and He + data is most sensitive to the assumed shock 
obliquity 6 Bn i (among other input parameters), the observational uncertainty 
in 0Bni maps over to an uncertainty of around ±1.5 in the value of i] . 

There are clear differences between the simulation results and the obser- 
vations. The He ++ distribution of thermal solar wind ions appears slightly 
narrower than the published observations. More notable though is the fact 
that the accelerated thermal He ++ ions are injected somewhat less efficiently 
in the simulation than in the observations, an inefficiency characteristic of 
highly oblique shocks that is present also for the solar wind protons, though 
not explicitly displayed in Fig. 1. The efficiency of acceleration of thermal ions 
could be increased via several means: (i) by lowering the shock obliquity angle 
9 Bnl , for which there is a large observational uncertainty; (ii) by decreasing 
rj , corresponding to increased turbulence, without altering the pick-up ion 
acceleration efficiencies substantially, and (iii) increasing the temperature of 
the thermal ions somewhat, though this would reduce the compression ra- 
tio and accordingly steepen the non-thermal continuum. However, such pa- 
rameter adjustments only incur small changes to the widths of the thermal 
distributions. Trial runs indicate that an agreeable fit can be obtained to 
the entire a particle distribution simply by reducing rj to rj ~ 1 , without 
changing other parameters. It is not clear what plasma characteristic would 
establish a species-dependent i] , a circumstance that contrasts the inference 
of r](H + ) « 7]{He ++ ) in the 91097 event analyzed by Baring et al. (1997). 

An instructive diagnostic on the acceleration model is to probe the spatial 
scale of diffusion. This is performed by examining upstream distributions of 
high energy particles. In their 1994 paper, Gloeckler et al. discussed an energy- 



9 




291 292 293 

Log 10 [Day of 1991] 



Fig. 2. The flux variations of accelerated pick-up ion populations as a function of 
time near the shock of Day 292, 1991. The data for 5 keV and 200 keV pick-up 
H + are depicted by filled triangles and squares, respectively, and are taken from 
Gloeckler et al. (1994). The Monte Carlo model generated fluxes at different dis- 
tances normal to the shock, and were converted to spacecraft times by incorporating 
solar wind convection. The 5 keV and 200 keV pick-up H + traces are displayed as 
dashed and solid curves, respectively, and exhibit an exponential decline upstream 
of the shock that is characteristic of diffusive shock acceleration. The model nor- 
malization protocol is discussed in the text. The lightweight line labelled "100 keV" 
is an approximate indication of a model prediction for the expected flux variation 
upstream for 100 keV protons. 

dependent rise in fluxes of non-thermal particles prior to the shock crossing. 
This was cited as indicating the existence of a pre-acceleration mechanism. 
The Monte Carlo simulation was run with flux measurement planes placed 
upstream of the shock at different distances, as well as downstream. This 
enabled determination of the number fluxes of a particular energy ion at dif- 
ferent positions upstream of the shock. This was converted to observer's times 
using the approximate speed of the Ulysses spacecraft in the solar wind. Re- 
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suits for two different H + ion energies are displayed in Fig. 2, together with 
corresponding data from Fig. 3 of Gloeckler et al. (1994) for identical energy 
windows. H + was chosen by Gloeckler et al. since good spatial resolution at 5 
keV was possible. Observe that the y-axis in Fig. 2 represents ratios of flux to 
time- averaged flux at the respective energies, so that the normalization does 
not explicitly exhibit the property that the absolute flux at 200 keV is always 
much less than that at 5 keV when downstream of the shock. Note also that 
the Ulysses data normalization was established by averaging over 3 days of 
accumulated data, whereas the model normalization was adjusted to match 
observed fluxes around 1/2 day downstream of the shock. 

The essential result of this comparison is that the spatial scale of the expo- 
nential decline of ions upstream of the shock is more or less identical to that 
of the model, for our choice of rj = 5.5 . Theoretically, high energy particles 
with a mean free path according to Eq. (1) establish an exponential dilution in 
space/time due to random scattering of the particles seeding upstream leak- 
age against a convective flow. From standard diffusion-convection theory, the 
spatial scale for this dilution is (k\\ cos 2 Obui + sin 2 0Bni)/ui.x f° r a shock 
of obliquity Qbu\ , where k\\ — Xv/3 is the component of the spatial diffusion 
coefficient along the mean field, k± = k\\/(1 + rf) , and A is prescribed by 
Eq. (1). Hence, this upstream dilution scale is proportional to the proton's 
energy. For the 200 keV ions, the simulation results are clearly well correlated 
with the data prior to the shock, modulo plasma fluctuations. On the other 
hand, for the lower energy 5keV ions, the exponential decay has a very short 
time scale, a factor of 40 smaller than for the 200 keV ions, and drops to back- 
ground levels very quickly. So, although the simulation results are consistent 
with the observed results, it is impossible to draw more definitive conclusions 
without an improvement in data time resolution. The 200 keV ions, however, 
have much longer mean free paths and travel much farther upstream, and are 
well correlated. 



4 Conclusions 

The comparison of phase space distributions from the Monte Carlo sim- 
ulation of diffusive shock acceleration with those observed by SWICS and 
HI-SCALE in the Day 292, 1991 event reveals a good deal of consistency be- 
tween theory and experiment for both the proton and He + spectra above 
speeds around 600 km/sec. At these speeds, the injection of pick-up protons 
dominates that of solar wind protons, and the acceleration of solar wind He ++ 
is inefficient relative to that of pick-up He + . The characteristic of the shock 
most relevant for these high energy components is its compression ratio r : to 
match the observed steep spectra, r ~ 2.1 is required, on the low end of the 
range quoted by Gloeckler et al. (1994). 
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The normalization of the energetic ion power-laws is best fit with a model 
"turbulence" parameter r] = 5.5 ± 1.5, where A = f]r g , corresponding to a 
ratio Kj_/ K\\ ~ 0.03 of components of the spatial diffusion tensor. For such 
values, cross-field diffusion is insufficient for the Monte Carlo model to account 
for observed accelerated He ++ , which requires rj ~ 1 so that the Bohm limit 
k± ~ k\\ is approximately realized. While the high energy protons and He + 
are modeled fairly well, there are significant discrepancies in the high velocity 
wings of the thermal protons. These deviations are slightly more marked than 
those found by Baring et al. (1997) in IP shocks without significant pick-up ion 
components. Clearly some non-diffusive element of heating may be present in 
the plasma shock that is not incorporated in the simulation. To explore such 
possibilities, it is planned to incorporate cross-shock and shock layer charge 
separation potentials in the Monte Carlo code in future extensions, to effect 
a more accurate modeling of the shock structure. It must be noted, though, 
that uncertainties in the shock parameters significantly impact inferences of 
discrepancies between theory and experiment. 

The flux increases of energetic protons seen upstream of the shock are quite 
consistent with the expected upstream leakage associated with diffusive shock 
acceleration. The value of rj = 5.5 inferred from the spectral fit scales the 
upstream diffusive lengthscale, and the accompanying exponential decline in 
predicted flux is commensurate with the Ulysses data presented in Gloeckler 
et al. (1994). Hence, the observed upstream flux precursor is not clear evidence 
of a pre-acceleration mechanism, as claimed by Gloeckler et al., though it is 
quite possible that some pre-acceleration mechanism may be acting. A most 
enticing feature of this work is that a single model parameter can couple 
injection efficiency to the spatial scales of ions upstream of the shock. Whether 
or not this identifies diffusion as the dominant operating mechanism remains 
to be seen, but it is clear that processes acting on the scales of a few gyroradii 
control both injection and upstream transport in this interplanetary shock. 
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